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Abstract 

We review Biichler and Colangelo's result that leading divergences at any loop 
order can be calculated using only one-loop calculations and we provide an alternative 
proof. We then use this method to calculate the leading divergences of and thus the 
leading logarithmic corrections to the meson mass in the massive 0{N) nonlinear 
sigma model to five-loop order. We also calculate the all-loop result to leading order 
in the large N expansion by showing that only cactus diagrams contribute and by 
summing these via a generalized gap equation. 



PACS: ll.lO.Hi, ll.15.Pg, ll.30.Qc, 12.39.Fe 



1 Introduction 



In a renormalizable theory the use of Renormahzation Group Equations (RGE) is common 
practice. RGEs have not yet received the same attention in non renormahzable effective 
theories. This is partially due to the fact that one does not normally have the problem of 
evolving coupling constants through a large energy range, and partially to the fact that 
RGEs in non renormalizable theories get more complicated as one goes to higher orders. 
Loop corrections however can be significant [U Ej [3] and need to be dealt with. 

Consider for example the scattering length Og in tttt S'-wave, 1 = scattering in ChPT. 
Close to threshold, this amplitude may be expressed in terms of the expansion parameter 
(Mn-/47rF7r)^ ~ 0.01. Despite the smallness of the expansion parameter, the one loop 
contributions cause a 28% corrections to the tree level prediction [5]. The reason for this 
is that beyond tree level the expressions for observables contain non-analytic functions [5] 
such as log M ^ which may be large if is small, even near threshold. It is only 
natural then to wonder about the size of higher order n contributions like (M^ logM^)", 
the so called leading logarithms, and about the size of their coefficients. 

In a renormalizable theory these coefficients are fully determined by a one-loop calcu- 
lation. This is a consequence of the RGEs The difference between a renormalizable and a 
non-renormalizable theory is that in the first case the counterterms needed at any given 
order have the same form, while in the latter case new ones are needed at every order. 
Nevertheless, one can still make predictions on the leading logarithms. This fact was first 
pointed out by Weinberg [1] in the context of Chiral Perturbation Theory (ChPT). He 
showed that at two loops the coefficient of leading logarithm (M^logM^)^ can be deter- 
mined simply by performing a one-loop calculation. 

This method has since been used in ChPT [U Ej |3] to two loop order in vrvr scattering 
[6] and in general [7j. It is no longer much used in the purely mesonic sector, since most 
processes are actually fully known at two-loop order as reviewed in e.g. [8]. 

In the last few years though, Weinberg's argument has received new attention [9l [TOl 
[TTl \T2\ [13], especially since Biichler and Colangelo were able to generalize the result to 
all orders [H]. They showed explicitly that one can obtain the (coefficient of the) leading 
logarithm at any order by simply performing one-loop calculations and that this coefficient 
is just a function of the lowest order coupling constants. The relevant part of their paper 
[9] and their algorithm to find the coefficients is described in Sect. |2l 

In Section 12.31 we provide an alternative proof of their results which does not rely 
explicitely on /3-functions but follows directly from the fact that all non-local divergences 
must cancel. This version of the proof has the benefit that it shows immediately that one 
only needs to calculate the divergent part without worrying about classifying higher order 
Lagrangians and that there is a direct link between divergences and leading logarithms. 

We then apply this method to the case of a massive 0{N) non-linear sigma model and 
calculate the corrections to the meson mass up to five loops. Sect. O contains a detailed 
explanation of our calculation. Similar calculations have been performed by [121 [13] who 
showed how to obtain the leading logarithms in the massless case by deriving a recursion 
relation for all possible vertices with up to four mesons. Since in the massless limit the 
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tadpoles vanish, this allows obtain the leading logarithms in a straightforward fashion. 

The authors of [TOj instead calculated the two-point function up to five loops in ChPT 
in the chiral limit using dispersive method^]. Once the first five leading logarithms where 
known, the next step was finding an algorithm that would allow them to calculate the 
n-th order one and eventually to resum the series. In paper [TTj they considered a linear 
sigma model and compared the correlator leading logarithms they found with those from 
ChPT, both in the chiral limit. They showed that it is not possible to simply use RGEs 
in the linear sigma model to resum the chiral logarithm series. The two scales present in 
the linear sigma model both generate logarithms that cannot be disentangled. 

We are however able to calculate the the five-loop meson mass in the massive non- 
linear case. What allows us to obtain these results is the observation that while one needs 
a complete Lagrangian to any order, this does not have to be minimal (see Sect. 12. 2p . nor 
does one need its explicit form in terms of chirally invariant operators. This, combined 
with the power of FORM [14J , allows us to calculate the relevant parts to five- loop order 
for the meson mass. 

A calculation of this magnitude needs as many checks as possible to ensure that no 
mistakes are made. We use two main checks. We perform the calculations in differ- 
ent parametrizations of the 0{N) nonlinear sigma model. Since different parametrizations 
distribute contributions very differently over the the various Feynman diagrams, the agree- 
ment provides a rather stringent check. Another check is that the leading term in N agrees 
with the result of the large N expansion of the model. 

To the best of our knowledge, this is the first time a study of the massive 0{N) is 
performed in this limit, whereas there is a vast literature on large in the massless case, 
especially for the linear sigma model, see [15| and papers citing it. 

We solve the mass in the 0{N) model in a fashion similar to that often used in the 
Nambu-Jona-Lasinio model, see [TB] and references therein. We first proof that only "cac- 
tus" diagrams contribute and that they can all be recursively generated by an equation for 
the exact propagator, a generalization of the usual gap equation. Finally, we are able to 
find a compact implicit expression for the all-loop physical mass. A thorough explanation 
of the study can be found in Section HI 

This paper is organized as follows. In Sect. [2] we discuss the results of P| (Sect. 12. ip 
and use their proof to show that one does not need a minimal nor a fully symmetrically 
formulated Lagrangian at higher orders (Sect. 12.20 . Sect. 12.31 provides an alternative proof 
of the results of [9]. Sect. [3] defines the 0{N) nonlinear sigma model and the different 
parametrizations that we use. Section S] discusses the large A^ case to all orders. The 
calculation of the leading divergences for the mass to five-loop order in the general case is 
described in Sect. |5l there we also give the result of the calculation. Sect. [6] summarizes 
our main results. The way we perform the integrals is described in App. |A1 

^In the same paper the authors also calculate the dispersive part of the three-loop pion form factor. 
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2 Renormalization group arguments 



2.1 The equations for the divergences 

This subsection recapitulates the parts of [9] we will use. We use dimensional regularization 
with d = A — w. The Lagrangians can be ordered in an expansion in h. We denote the 
lowest order Lagrangian with Cq. 

n>0 

= + (1) 

The divergent part contains the inverse powers of w needed for the subtraction of the loop 
divergences at order n [§. 

We now expand each term into a set of A''„ operators ol""^: 



i 

i=l 

N„ / n An) 



= ElE^l^l"^- (2) 



i=i \k=i 



One key difference between a renormalizable Lagrangian and these is that here the 
minimal basis of operators Of^^ grows with n. In the remainder, we will assume that 
all one particle irreducible diagrams (IPI) are made finite separately. This simplifies the 
calculations and arguments. That this can always be done is discussed in [9]. We have 
already used here, in the expression of £^'^, the fact that all divergences are local. The 
Cj-""^ are usually referred to as Low- Energy-Constants (LECs). 

The bare Lagrangian is fi independent. This leads to the equations 

= /i-^/:^"--^ 
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(3) 



These must be fulfilled separately for each ra, i and inverse power of w. We define the /3 
functions via 

d (n) n{n) (n) , , \ 

The RGEs of ([3]) thus become 

n / rl \ A^"^^ 

A'"'.g(-.„. + ,A):^.o. (5) 



^ When comparing with ChPT one should remember that our order ?i" corresponds to the order p2n+2 
in ChPT. 
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The A^^^ do not depend explicitely on the fi. They do, however, depend upon the c\^\ since 
they must cancel the divergences stemming from the loops of the lower order Lagrangians. 
We can thus simplify ([5]) to 



PI 



(n) 



E 

k=l 



-nw + ^ 



( im) 

mwCj 



d 



(m) 



0. 



(6) 



The loop contributions must be polynomials in the coupling constants c\"^^ . It follows that 
A^^^ I must also be a polynomial. This means that we can split further into separate 
equations. It must be true for each power in w, but also for each c\^^ monomial, since 
these are independent parameters and can, in principle, be varied freely. The coefficient 
of each separate monomial in the c-"^ in then must vanish separately. It is this extra 
information that allows us to obtain all the leading divergences from one-loop calculations 

The powers of h in any diagram come from two places, a factor comes from the 
number of loops / and the remainder from the h"'Cj^^ present in the diagrams when vertices 
with ol""^ occur. This shows that each monomial in the c-"^ will also come from a well 
specified loop level. 

Let us write out the various equations for the first few orders. To order h, we only get 
the equations 



PI 



(1) 



0. 



(7) 



There is no dependence on any of the higher order LECs. We now introduce a notation for 
the LEG dependence. For both /^^"^ and A^^i^ we add a subscript / indicating the loop level 
it came from and an argument indicating its polynomial dependence on the Cj^^. In order 
to simphfy notation in the following we shall omit the i index, so that now /^f"^ and A^i^\j) 
become p\^\j) and A^^ , but one should remember throughout the calculation that we are 
always speaking about the component /3f or v4[^j- . At one loop this only adds a subscript 
1 



P^' 



At order the /3| function and can have a first order dependence on c^-^ , we make 
this explicit 



aS^ 

Af^ 



(9) 



A sum over ji is implied. The ji indicates that one should consider the minimal set of 
operators available at order 1. Putting ([9]) in ([6]) gives three conditions on the various 



4 



coefficients of the monomials in the c-"^ 

/^f = 24?, 

24? = (10) 

At order the polynomial dependence gets more complicated 

A?^ - ^g^+4;Ug)(jO + cgUff(2j) + c«4US?(j,A:0. (11) 

Which terms can show up at which level follows from the h counting and the fact that / 
loops can at most diverge like 1/wK Putting fill I) in gives the relations from the 0{w^) 
in dSD 

Pi'' = 34?, 
/3f(ji) = 24?(ji), 

P?\j2) = 4¥(J2), 

P?\j,h) = 4?(jifci)- (12) 

The —nw and J2m,j''^'^'^f^'9/dc^j^^ terms in (jSj) always combine to give exactly the loop 
level back, see [S] for the general proof. This gives the first equation in f[T7l) . 
The part gives 

34? = /5S4?(ji) + /5g4?(j2), 
242(j\) = 2/?;i;4?(jifco + /3a(ji)4?(A:2). (13) 

Here we have written out the operator C-"'' subscripts (i,j2, ^i, ^2)- That is to stress the 
fact that there is a sum over a different index ji 7^ i. We will omit it in the following. 
In deriving the second equation we have used that ^nHJifci) is symmetric in jiki and 
relabeled some indices. The final equation comes from the w~'^ part and reads 

34?=/5S4?(ji). (14) 

The set of equations gives 

64? = [2/5;i;^s?(j\A:o + pfiijMfhh)] . (is) 

As one can see, the leading divergence can be calculated with purely one-loop calcu- 
lations. 
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The argument above can be generalized to all orders, see [9], here we only quote the 
results. We introduce the notation 



v. = EE/^rT^ (16) 



m=l J UCj 



The general set of equations thus reads 

f^i '* ~ ^A\i ^ ) / = 1, . . . , n , 

l~k+l 

l^t^ = E '^i'Ai-i',k-i , l = k,...,n; k = 2,...,n. (17) 
i'=i 

For the leading divergence A^^n this equation reads 

= 4?- (18) 

which is a generalization of (171), ( ITOl) and ( ITSl) . This is the main result of [9j we will be 
using. This is a recursive relation. The order n = 1 counterterm A^iiO'^^^ is fixed by the 
requirement that it should cancel the Cq one loop 1/w pole. This is now the order n = 1 

(2) 

coupling in the Lagrangian. The order n = 2 coupling A22 is fixed by the requirement that 
it cancels the one loop l/w"^ pole coming form the £1 with the y4.ii'' coupling we fixed in the 
previous step. And so on, the n-th order A^J^ is fixed by the requirement that it cancels the 
1/w"' divergences. These can be calculated by considering all one loop diagrams generated 



n. 



by the A^. _i „_iC^""^) ,. • • , -AiVO(^) and £0 vertices that can contribute to order 

The renormalized coupling c^^^ now contains a log/i. This exactly cancels the log/x 
dependence that comes from the Cq loop integral, which has the same coefficient as the 
divergence. Analogously, when an observable is calculated up to order n its expression 
contains a (log/x)" term whose coefficient is given by the A^^l So (HBD gives a recursive 
expression for the coefficients of the leading logarithms. 



2.2 Nonminimal sets of operators 

In Section [2?T] we shortly went through the arguments of [9j to derive the leading divergence 
at any order from only one-loop calculations. In practical applications of the formulas 
above one needs a classification of the terms O^""^ needed at each order n. Determining the 
complete and minimal set is in general rather complicated, see e.g. [17J for C2 in ChPT. 
Luckily we do not have to have a minimal and complete Lagrangian in general. 

It is sufficient to have a Lagrangian that is complete for the particular process at hand 
and lower order Lagrangians that are complete enough so that all needed /9f^(l) can be 
obtained. The Lagrangian does not need to be minimal since the arguments in Section 
12.11 relied on the fact that all c\"'^ can be varied independently. If we add an irrelevant 
term, e.g. one that vanishes via partial integration or other identities, its coefficient can 
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definitely be freely varied and will not show up in any actual higher order calculations. 
This also means that the f3 function of this irrelevant term can be chosen freely since it 
will never appear in any expressions. 

If we now have two related terms c^"''(9^"'' and c^^^O^\ we can always write them as an 
irrelevant one cj^^Oj-"^^ and a relevant one c^J^jol^j. Since the irrelevant combination will 
never appear and its f^lj^J function is free, we can just as well leave both terms Ci^^Oi^^ 
c^2^02^^ in the Lagrangian and leave the (5^^^ spread over the original two terms Z^^"'' and 
/32"^ How it is split between the two terms depends upon the choice of A^^'J'', which is free. 

In practice, it is sufficient to calculate the divergences and to express them in some 
complete set of operators As long as the set is complete for the given application we 
will obtain the correct result. 

When constructing a Lagrangian one normally takes into account all terms that have 
the correct symmetry at the required order and then one removes the so called equation of 
motion terms, see [I7j App. A, for a discussion. In this case, however, we will use neither 
constraint. We will keep the equation of motion terms, since then we can make all IPI 
diagrams finite, see the discussion in We will also use a standard Feynman diagram 
calculation to obtain the infinities and not a more sophisticated method such as e.g. the 
heat kernel expansion that was used in [T8] . 

The reason is that this way we can use standard Feynman integral techniques and we 
do not have to evaluate all the divergent combinations of propagators that can appear, see 
[T8| [19] . which is rather difficult at higher orders. The drawback is that this procedure 
breaks the symmetries of the Lagrangian in individual parts of the calculation even though 
the final result will respect all symmetry properties when we use dimensional regularization. 
However, as we saw in the previous section, we only use the divergent parts of these terms 
and these must obey the symmetries since they are recursively determined by a symmetric 
lowest order Lagrangian. 

The answer for the divergences for a given process will thus be correct even without 
explicitly fixing counterterms with Ward identities. The correct combinations must show 
up in our procedure. 

The conclusion from this section is that we simply calculate all one-loop diagrams and 
rewrite them as terms in the Lagrangians, without bothering to check if we have a minimal 
Lagrangian. 

We also do not need to have a complete Lagrangian, an operator O^^'' will only be 
relevant if it has a nonzero (3 function. We thus let the calculation itself produce all terms 
that have a divergence, give them a coefficient cf'\ and use those in the equations derived 
in Section 

For the subleading divergence, the same type of argument shows that it is sufficient 
to have the lowest and first order Lagrangian in a symmetric form to get the subleading 
divergence at all orders and the obvious generalization to the further divergences. 
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2.3 An alternative proof 

A more direct proof of the results of Biichler and Colangelo is also possible. 

We have presented their method as well since the arguments in Sect. 12.21 made use of 
their formulation of the proof in [9j. 

We rely here on using only IPI diagrams and assume they are made fully finite, as 
was shown to be possible in ^ . In this section we obtain the same relations in a more 
transparent fashion. We first rewrite ([1]) and ([2]) as 

c"^- = y: hy— E f E ol-^ . (19) 

?i>0 i \k=0 / 

We introduce the notation {c}" to indicate all possible combinations c^k^j^c^j^j^ ■ ■ ■ ci™v '^i^h 

TTii > 1, such that J2i=i,r^i = ^ ^-iid J2i=i,T = I- The c^"'' with k >1 have no direct /i 
dependence. They only depend on ^ through their dependence on lower order parameters. 
The Cq^'' do depend directly on /x. Note that since we consider only IPI diagrams we have 
that {cYn = {ci?}. 

We denote the contribution from all /-loop diagrams at order as and we expand 
this as 

= E Llw-' . (20) 

fc=0 

This only includes the divergences coming from the loop integrations, not those from the 
coefficients in the Lagrangian. 

A main observation [9] is that a given loop level at a given order ff' always comes with 
the same power of /i because of the way the powercounting works. 

We can now study the contributions at the different orders in h and 1/w. For clarity 
we add here as well which combinations of couplings of order n > 1 the results depend on. 

At order ff' we have only Lq. At order we have 

^ [li--Ll,{{c}\) + L\,) + /."-4o({c}^) + L\, . (21) 



w 

The divergence must cancel so to get the divergent combinations we have that 

LU{c]\) = -L\,. (22) 

This allows to determine the divergences that need to be subtracted from a one-loop 
calculation and it shows that by expanding and taking it; — > the explicit log/i 
dependence of any process is 

-log/xL^o({c}l) =log/iLji. (23) 
At order the full contribution is 

+- [^^-''"LU{c]\) + ^l~-LU{c]l) + ^.-L\,{{c]\) + lI,) 

w ^ ' 



f28) 



All divergences must cancel, also those with powers of log /i. If we only look at the parts 
with l/w"^ and log/i/w we obtain two equations 

Ll{{c]l) + LU{c]\) + Ll, = Q, 

2LU{c}l) + LU{c}\) = Q (25) 

The difference in the coefficients from the first to the second equation comes from the 
expansion of the different powers of . These equations have the solution 

-^Oo({c}2) = -^22 5 

Ll,{{c}\) = -2Ll,. (26) 
The leading logarithm can be obtained by expanding in fl2^ and using fl26|) : 

llogV(4Lgo({c}^)^' + ^n({c}i)) =logV^^2- (27) 

So here we reproduce the known result and that it can be obtained from a one-loop cal- 
culation. The calculation for the next two orders follows the same lines. A clear pattern 
emerges. 

At order the leading part is given by 

^ (/i-""^00({c}::) + LUicYnZl) + ■■■+ /i-'^^^l n-l{{c]\) + 

This part is the one that contributes to the logyu/iy"'~\ . . . , log"~^ ^/w divergences 

leading to the set of equations: 

n 

^^^L^__,({c}i)=0 j = 0,..,n-l. (29) 

with 0° = 1 and L"„({c}o) = The generalization of the solution then is 

^L.n-.({c}D = (-ir(^)^L- (30) 

We can prove that this solves the equations (l29l) above by observing that they can be 
written as 

Hi3(|:(4)'''""^s(w-') = °- 

Plugging in (130|) we see that this becomes 

lim (^a^y (-a + 1)"L1 = , (32) 

which is clearly satisfied. Using lim^^i {a-^ (—a + 1)" = [—l)^n\ one can also derive 
that the dependence on log"/! is 

log>Ll. (33) 
This completes our alternative proof of the main result of [9]. 
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3 The 0{N) nonlinear sigma model 



The 0{N + 1)/0{N) nonlinear sigma model has as Lagrangian 

F2 



(34) 



$ is a real + 1 vector that transforms as the fundamental representation of 0{N + 1) and 
satisfies the constraint <l>"^$ = 1. The second term is the one that breaks the symmetry 
explicitly by setting 







The vacuum is given by 



($^) = (1 0...0) 



(35) 



(36) 



which breaks the 0{N + 1) spontaneously to 0{N). There is both a spontaneous symmetry 
breaking triggered by the vacuum ( l36i) and an explicit one given by 

This Lagrangian corresponds to the lowest order Lagrangian of two- (quark-) flavour 
Chiral Perturbation Theory for = 3 [201 12] and has been used to describe alternative 
Higgs sectors in several beyond the Standard Model scenarios. 

As mentioned in the introduction, we make use of different parametrizations to check 
the validity of our results. We write $ in terms of a real iV-component |f| vector 0, which 
transforms linearly under the unbroken part of the symmetry group, 0{N). We use here 
four different ways to do this parametrization 



$1 



$9 




(37) 

(38) 
(39) 

(40) 



$1 is the parametrization used in [2], $2 the one originally introduced by Weinberg 
$3 is such that the explicit symmetry breaking term in flMl) only gives a mass term to the 
field but no vertices. $4 is the parametrization one ends up with if using the general 



^We refer to these as a flavour components. 
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prescription of [2T|. These are all examples of the parametrization that keeps the 0{N) 
symmetry manifest: 

Here f{x) is any function with /(O) = 1. 

One thing is worth mentioning. In this work we always calculate with the usual Feyn- 
man diagram techniques. In our calculation we split the A^-vector field into an external 
{(^e) and a loop (^) field, 4> ^ 4>e + C The divergence structure we obtain is expressed in 
terms of (pE and we then set 4>e ^ 4> and use that as input for the next step. Splitting 
in this way, the symmetry is no longer manifest in each term. Since we are renormalizing 
these terms, it would be nice if they were obviously symmetric. We could have used the 
background field method, used e.g. in [21 [18], and split into a classical and a quantum 
field with well defined symmetry properties and then calculated the divergent part up to 
a given number of external legs. This way we would have been assured that our divergent 
Lagrangian can be rewritten into terms fully obeying the symmetry. Once that is done, we 
could then use the symmetric quantum field again for the next step. The problem is that 
the rewriting into symmetric terms is not easy to implement. We have shown in Sec. 12.21 
that our method gives the correct answer too. 

One additional way to check that our results are correct is to compare with the known 
results. For = 3 our 0{N) corresponds to SU{2) x SU{2)/ SU{2), and loop corrections 
to the pion mass up to two loops are fully known in this model [221 [23]. The leading terms 
at two-loop order were first obtained in [6j. They read 

Lm = T^log;^^- (42) 
where in numerical applications one usually chooses Ai = M. 



4 The large N approximation 

The linear sigma model has been treated very much in the large approximations. The 
literature can be traced back to \i5\ but the literature on the nonlinear sigma model 
is smaller. In addition, it is mainly restricted to the massless case while here we are 
interested in the massive case. [21] did include masses but only to first order. There are 
some subtleties involved in large N in effective theories because of the presence of the 
higher order Lagrangians, see e.g. [25]. We however keep our discussion on the level of 
the loop diagrams with the lowest order Lagrangian and stick to in the notation of 
Sect. [231 

We choose here the Lagrangian to be extensive in A^. This means we have to choose 
oc A^. Looking at the Lagrangians in Sect. [3] we easily see that vertices with 2n legs 
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Figure 1: A typical diagram that contributes at leading order in N. Note that vertices 
can have many different loops attached since the Lagrangians contain vertices with any 
number of fields. The flavour-loops coincide with the loops in momentum. 

have a factor _p2-2n ^^^^ thus suppressed by A^^"". Extra factors of come from closed 
loops in the flavour index. It is sufficient to look at one-particle-irreducible (IPI) diagrams, 
the flavour indices for lines not inside a loop are determined by the external flavour indices. 

We thus look at the one-particle-irreducible diagrams only and use methods similar 
to those used in [1] for proving the powercounting. A given diagram has Ni loops, 
vertices with 2n legs, Nj (internal) propagators and Ne external legs. These are related 
via 

n 

2Nj + Ne = Y.'^nN^n. (43) 

n 

We thus get 

Nl = - l)^2n -\ne + 1. (44) 

The tree level one-particle-irreducible diagram has one vertex with 2n = Ne and thus 
comes with a factor of N'^-^^/^. A generic one-particle-irreducible loop diagram has thus 
a suppression factor 

from the vertices. Extra factors of A^ come from the closed flavour loops where each closed 
flavour loop give a factor of A^. (145!) shows that for a diagram to be leading in A^, there must 
be as many closed flavour loops as there are loops. Since the fields only carry one flavour 
index, this means that each loop must coincide with the flavour loop and there can be no 
lines shared between two loops. All diagrams that contribute to a given process at leading 
order in A^ are thus diagrams that only contain products of one-loop diagrams, these we 
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Figure 2: The graphical representation of the equation that generates all the cactus dia- 
grams for the propagator. A thick line indicates the full propagator, a thin line the inverse 
one. 

call cactus diagrams after the looks of the "prickly pear" cactus. A typical example is 
shown in Fig. [H 

How can we now resum all of these diagrams? The easiest way is to notice that they 
can be generated recursively. First we note that the inverse of the full propagator is given 
by the inverse of the lowest order propagator and the sum of all the one-particle-irreducible 
diagrams with two external legs. This leads to the equation graphically depicted in Fig. [2l 
The difference with the usual gap equation in NJL-type theories as discussed in e.g. [16] is 
that we have an infinite number of terms here. This actually turns out to be manageable. 

Let us look at the Lagrangians of Sect. [3l They all contain at most two derivatives. 
The loops in Fig. [2] are all tadpoles and thus produce no extra dependence on the external 
momentum p. All the dependence on p must come from the derivatives present in the 
vertices. The full inverse propagator is thus of the form Zpp'^ — Bp, where neither Zp nor 
Bp depend on the momentum p. This is true for all parametrizations. 

We now use the first parametrization. The vertices with derivatives are generated by 

^ ^ r0'a^0V'^5^</'% (46) 



2FH- 



F2 



where we have written (f)^(f) = to bring out the sum over flavour indices explicitly. Each 
loop must allow for a sum over the flavour indices to be leading in A^. The derivatives 
must either both act on the external fields or both on the same loop to give a nonzero 
resulll^. When they act inside a loop, the fields d^cjy^d'^cl)'^ must be contracted to form the 
loop. This requires h = c and the flavour in this loop is thus determined by the outer fields 
and cannot be separately summed over. Consequently, the diagram gives no leading A^ 
contribution. 

If the derivatives hit the external legs, the indices in at least one loop are fixed by the 
external ones and again cannot be summed over. It follows that the contribution is not 
leading order in A^. 

In either case the loop diagrams generated by the kinetic term give no leading A^ 
correction, so = 1. 



This will not be true for more complicated processes but can be dealt with in that case as well [26] 
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Thus we only need to look at the vertices coming from the mass term 



C^ass = F^My 1 - rr/F^ = F'M'fix) = F'M' <^^^' , (47) 

i 

with X = (f)"'(f)^/F'^. Again, consider the a loop diagram. The external legs need to come 
from the same flavour index otherwise it will not be leading in For each term in (j4j 
there are i ways to choose which x corresponds to the external legs 



i>2 



The sum in (1481) starts from i = 2 since the i = 0,1 terms are a constant and the tree level 
mass term respectively. Eq (H7I) then reads 

-^M20L0LJ^^=^-ll • (49) 



F2 



The leading contribution comes from contracting the fields with the same flavour index. 
There is only one way to do this for each term. Each contraction corresponds to a tadpole 
in Fig. [21 The full result can be written as 



- Bp = - - \ I ^ _ = - 1 I , (50) 

F2- 



^1 + 4^(5, 



where iA{Bp) = J d'^pl/{p^ — Bp) is the relevant one-loop tadpole integral. The all-loop 
result at leading order in is thus the solution of 



= M^^^ 1 + ^AiM^s) ■ (51) 

Here Bp coincides with the physical mass squared, M^f^y^, since Zp = 1. 

The same result can be derived in the other parametrizations. If we take the third one, 
where the only vertices come from the term with derivatives, the same type of argument 
as above with derivatives and flavour indices shows that the only relevant vertex is 

^ (f)ydf,(f)^dy'' . (52) 



8F2 



So here the gap equation reduces to the first nontrivial term on the right-hand-side only. 
The structure of the inverse propagator is still Zpp"^ — Bp and the gap equation leads to 
two equations with M^f^ = Bp/Zp 



1 + J^MM; 



2 



AF'^Zp^ ^''y ' 



ZpMl^y^ = - ^^M^^^AiM^J . (53) 
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Here we have expressed the integral containing an extra using (IMI) in terms of the one 
without. Solving leads to the solutions 



1 / ^ ^ 

2 \ ^phys, 

N 



M' = M^,^, (l + j^AiM^s)) , (54) 
which agrees with the previous result (j5Ti) . 



The leading logarithm can be expressed by replacing A{Mp,^y^) by 
In terms of 

we can invert the result (|5T1): 



M2. 1 So 3 231 4 7 . 7293 ^ . 1062347 o 

^ = l--,H--,= -,» + _,*--/ + _/-15y'H-^^y» + .... (57) 

Note that (15T|) actually converges faster. Expanding the square root in 

,-^!^§-,0K^ (58) 

we have 

M2 ^1 1 2 1 s 5 4 7 . 21 33 7 429 o 



Mp\^, 2 8 16 128 256 1024 2048 32768' 

which has much smaller coefficients than fl57p . 



5 The calculation 

We determine the mass by finding the zero of the inverse propagator. We therefore must 
calculate all the one-loop diagrams that are needed to obtain the divergence of the inverse 
propagator to the order desired. At order h there is only one diagram, at h"^ there are 2 
and at order there are 4. These are shown in Fig. [31 We have not shown them but at 
order h'^ there are 7 and at there are 13 diagrams to be calculated. 

To order h it is sufficient to know the lowest-order Lagrangian, but at order h we need 
to know the (divergent part of the) vertices coming from the Lagrangian of order h with 
two and four external legs. The diagram of Fig. [3]^a) gives the divergence of the vertex 



15 



I 
la) 




Figure 3: The diagrams needed up to order 3 for the inverse propagator. Vertices of order 
are indicated with |T]. (a) The diagram needed at order h. (b) The 2 diagrams needed 
at order h'^. (c) The 4 diagrams needed at order h^. 





Figure 4: The diagrams needed for the divergence of the 4-meson vector, 
diagrams to order h. (b) the 4 diagrams to order fi^ 



(a) The 2 



with two legs but we also need to calculate the divergence of the vertex with four legs. 
This requires the diagrams shown in Fig. Hl^a). 

To order h^, we need still more vertices, we need the divergence of the two-leg vertex 
to order fi^, these diagrams we already have but we also need the four- leg vertex to order 
h'^ which can be calculated from the diagrams in Fig. Il](b). Inspection of the vertices there 
shows we already have all we need but for the 6-leg vertex at order h. To obtain that we 
also need to evaluate all diagrams shown in Fig. [51 By now, the pattern should be clear. 



to get the mass at order we need the 2 and four- meson vertex at order , the 2, 4 
and 6-meson vertex at order h^~'^ and so on. Similarly one can see that to get the mass 
at order h"', we need to calculate one- loop diagrams with up to n vertices. The extension 
to order shows that we we need to calculate 18, 26, 33, 26 and 13 at orders h^, . . . 
respectively. 

We have organized this calculation by first expanding the lowest-order Lagrangian to 
the order needed, up to vertices with 12 mesons for this work. With these vertices we then 




Figure 5: The 3 diagrams needed for the divergence of the 6-meson vector to order h. 
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i 


tti for = 3 


tti for general 


1 


-1/2 


1 - 1/2 


2 


17/8 


7/4 - 7/4 + 5/8 A^^ 


3 


-103/24 


37/12 - 113/24 A^ + 15/4 A^^ _ jys 


4 


24367/1152 


839/144 - 1601/144 A^ + 695/48 A^^ - 135/16 A^^ + 231/128 A^^ 


5 


-8821/144 


33661/2400 - 1151407/43200 A^ + 197587/4320 A^^ 
-12709/300 + 6271/320 A^^ - 7/2 A^^ 



Table 1: The coefficients of the leading logarithm LJ^^ up to z = 5 for the physical meson 
mass. 

calculate all IPI diagrams with up to 10 external legs. The divergent part of all needed 
integrals can be calculated relatively easily using the technique described in App.jAj At this 
stage, the dependence on external momenta is also put back as derivatives on the external 
legs and everything assembled to give the divergent part at order h for all the vertices with 
up to 10 legs using fl22|) . So we have assembled everything we need to calculate the one- loop 
divergences to order . The 26 diagrams are evaluated and we obtain the divergences at 
order using fl26l) . The process is then repeated up to order . All of the above steps 
have been programmed in FORM. The CPU time needed increases rapidly with the order 
n one wishes to reach. The Lagrangians at higher orders tend to contain very many terms 
and constructing the diagrams with many external legs at higher orders is also extremely 
time consuming. The CPU time used on a typical PC for the mass- divergence to order 
was approximately 0.1 seconds for 0.3 seconds for 11 seconds for , 700 seconds for 
ti^ and 30000 seconds for fi^ . These running times were achieved after several optimizations 
in the choice of routing the external momenta through the Feynman diagrams. 

We have performed the calculation for each of the four parametrizations shown in 
Sect. [3l Since each parametrization distributes contributions rather differently over the 
different Feynman diagrams, this provides a strong check on the consistency of the final 
result. 

We can from these divergences then obtain the leading logarithm. This leads to the 
result for the physical mass 

Khys = M^(l + a,LM + a2Llj + a^Ll, + a^Ll, + a^Ll, + ■■■), (60) 

where Lm is defined in (l42l) . The coefficients oi, . . . , as are give in Tab. [T]for A^ = 3 and 
general N. The result for A^ = 3 agrees with the known results for oi and 02 given in (H2|) 
and the leading term in A^ at each order agrees with the expansion of the all-order result 
in the large A^ expansion (1571) . Note that the large A^ approximation is definitely not a 
good approximation to the A^ = 3 coefficients. 

The result (1601) can be inverted and we find again a better converging expansion. 

= M^s (1 + + hLl,^^^^ + hLl,^^^^ + h,L\,^^^^ + hLl,^^^^ + ■■■). (61) 



17 



i 


bi ioT N = 3 


bi for general 


1 


1/2 


-1 + 1/2 N 


2 


-13/8 


1/4- 1/4 iV- 1/8 


3 


-19/48 


2/3 - 11/12 N + 1/16 


4 


-5773/1152 


-8/9 + 107/144 iV - 1/6 iV^ _ x/^g jys _ ^^i2s ^4 


5 


-3343/768 


-18383/7200 + 130807/43200 N - 2771/2160 A^^ _ 527/I6OO A^^ 
+23/640 iV4 + 7/256 A^^ 



Table 2: The coefficients bi of the leading logarithm L\,j^^^^ up to i = 5 for the lowest order 
meson mass in terms of the physical mass. 

The coefficients 61, . . . , 65 are give in Tab. [2] for A^ = 3 and general N. Just as the coefficients 
in (159|) are much smaller than in (1571) we see that the bi are much smaller than the Oj. 

In order to get a feeling of the size of these corrections and of the convergence of the 
series for the very relevant case A^ = 3 we have plotted them in Fig. [6l On the left side 
we see the result and on the right side the result for a value of F = 0.090 GeV 
amd fi = 1 GeV. Both clearly converge in the region shown and the inverse one clearly 
converges faster. 

6 Conclusions and discussion 

In this paper we have obtained three main results. 

First we provided an alternative proof for the results of [9] for the leading divergence 
at any loop order. Our proof relies of course on the same physical principles as the one in 
[9] but is simpler algebraically. 

Our second result is the large A^ expansion of the massive nonlinear 0{N) sigma model, 
we did not use the analog of the methods in [15J for the massless case but obtained a recur- 
sive equation to sum all the relevant cactus diagrams. This method is clearly extendable 
to other processes than the mass we considered here. The formula obtained for the mass 
in the large A^ limit (15T]) is remarkably simple. Its relevance for the real case A^ = 3 is 
not clear, since the large A^ result gives a rather poor approximation to the Oj coefficients 
when A^ = 3, as shown in Tab. [H 

The third result is the actual calculation for general A^ of the leading logarithm for 
the meson mass to five-loop order. This result agrees with the known two-loop result for 
A^ = 3 and with the derived result for the leading term in A^. 

It is clear that the methods developed in this paper can be applied to other processes as 
well, both the large A^ method and the leading logarithms to higher loop orders for general 
A^. Work is in progress for the decay constant, vrvr-scattering and the formfactors [26] . 
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1.1 




0.9 I 1 1 1 1 ' 1 I ' 1 1 1 1 1 

0.05 0.1 0.15 0.2 0.25 0.3 0.05 0.1 0.15 0.2 0.25 0.3 

[GeV^] Ml^y, [GeV^] 

(a) (b) 

Figure 6: The expansions of the leading logarithms order by order for F = 0.090 GeV, 
/i = 1 GeV and = 3. (a) Mp\y, in terms of Eq. ([60]) (b) in terms of M^^^y^ 
Eq. dSU). 
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A Integrals 

We can get rather high powers of momenta in our the integrals. Let us first look at integrals 
without external momenta. These are of the form 

'"^ I J {2t:Y (q' - M^T ■ 
This vanishes for m odd, and for even m we can us^ 



^ 2W2 (l + i)! ^ 



G^/il.../im 5'a»1A12 • • • 9fJ.m-ll-Lm ~^ 9pLlpL3 .... (63) 



^ The formula is only valid for d = A which is sufficient for our purpose. We derived it using recursive 
methods but it is probably well known in the higher loop integral community. 
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Where g^^, is the metric tensor. The right hand side of the last term consists of all possible 
ways to put the Lorentz indices on the metric tensor g^^fj,^ and is symmetric under all 
interchanges of the indices. After that we use recursively 




1 + 



'2-M2 



(64) 



to obtain terms either without propagators or without powers of q^. Of the resulting 
integrals only two have a divergent part 



with (i = 4 - 2e. 

We need to do know the divergent parts of one-loop integrals with up to 5 propagators 
for this calculation. This we do by combining propagators using Feynman parameters and 
then shifting the momentum variable to obtain integrals of the type flS21) . The Feynman 
parameter integrals needed are always simple polynomial ones. 

The above procedure can be programmed in FORM to work recursively. 
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